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Abstract 

We investigate the Manakov model or, more generally, the vector nonlinear 
Schrodinger equation on the half-line. Using a Backlund transformation method, 
two classes of integrable boundary conditions are derived: mixed Neumann/Diriclilet 
and Robin boundary conditions. Integr ability is shown by constructing a generat- 
ing function for the conserved quantities. We apply a nonlinear mirror image tech- 
nique to construct the inverse scattering method with these boundary conditions. 
The important feature in the reconstruction formula for the fields is the symme- 
try property of the scattering data emerging from the presence of the boundary. 
Particular attention is paid to the discrete spectrum. An interesting phenomenon 
of transmission between the components of a vector soliton interacting with the 
boundary is demonstrated. This is specific to the vector nature of the model and 
is absent in the scalar case. For one-soliton solutions, we show that the boundary 
can be used to make certain components of the incoming soliton vanishingly small. 
This is reminiscent of the phenomenon of light polarization by reflection. 
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1 Introduction 



In this paper, we consider the vector nonhnear Schrodinger equation (VNLS) on the 
half-hne. Historically, the model on the line for a vector field with two components was 
introduced and solved by the inverse scattering method by Manakov yj. The VNLS 
on the half-line can be seen as a generalization of the famous (scalar) one-dimensional 
nonlinear Schrodinger equation (NLS) along two directions: adding internal degrees of 
freedom and adding a boundary. It is worth mentioning that each direction is not just 
a trivial generalization exercise. Indeed, such considerations have been at the origin of 
major breakthroughs in the history of integrable systems. For instance, when the origi- 
nal Lieb-Liniger model [2j for spinless particles with contact interactions was generalized 
by Yang |3] to particles with spin, this gave rise to the now famous rational solution 
of the Yang-Baxter equation. Similarly, when studying integrable systems with bound- 
aries, Sklyanin unraveled new algebraic structures controlled by the so-called refiection 
equation [4j. Concerning the NLS equation, these questions have been studied rather 
well at the quantum level. Vector generalizations on the line [51 E] and on the half-line 
[7J have been investigated essentially from an algebraic point of view. At the classical 
level, a very general framework for dealing with NLS equations related to symmetric 
spaces has been presented in [8J. Although the VNLS equation is only a particular 
case of this general theory, it is still the object of active research: for instance, the full 
understanding of the factorization property of vector soliton collisions has only been 
achieved relatively recently (see for instance [H] and references therein). At the classical 
level again, the general question of generating boundary conditions of arbitrary order 
compatible with a general multicomponent nonlinear Schrodinger equation (expressed 
in terms of Jordan triple algebras) has been treated in [lOj. However, to the best of 
the authors' knowledge (and surprise), the problem of formulating the inverse scattering 
method with integrable boundary conditions to tackle for instance the question of inter- 
actions of vector solitons with a boundary has not been addressed yet. This is in sharp 
contrast with the scalar case where a rich literature has appeared over the years since 
the early attempts of adapting the inverse scattering method to initial-boundary value 
problems iJlj up to the very general method developed by Fokas (see e.g. the textbook 
|12|). in passing by the Backlund transformation method found in [T3| [TU [T5]. Recently, 
another idea [16] was revived in |17) and used to investigate soliton interactions in the 
presence of an integrable boundary by means of a (nonlinear) mirror image technique. 
As far as integrable boundary conditions are concerned, all these techniques are related 
and it seems that the central object is a special Backlund matrix that allows for the 
implementation of the mirror image technique. In the context of Fokas' method, this 
matrix evaluated at the location of the boundary is the matrix involved in the so-called 
linearizable boundary conditions. 

The goal of this paper is to introduce and study the VNLS with an integrable bound- 
ary in order to lay the foundations for a deeper understanding of integrable vector field 
interactions with a boundary. Section [2] contains known results for the inverse scattering 
method of the VNLS on the line. This allows us to introduce our notations as well as 
to collect useful results for the rest of the paper. Section [3] represents the core of this 
paper. We first derive two classes of integrable boundary conditions using the Back- 
lund matrix approach. Integrability is established by constructing an explicit generating 
function for the conserved quantities for the VNLS on the half-line. This is then used 
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to implement the mirror image technique that allows us to use the inverse scattering 
method results on the line to obtain the reconstruction formulae for the fields on the 
half-line. In section 21 we use our general results to investigate the behaviour of one 
soliton bouncing off a boundary. An interesting phenomenon of transmission between 
different modes of polarization is demonstrated. It is due to the interaction between the 
soliton polarization and what we can call the natural boundary polarization. This will 
be made more precise in the text. Our conclusion are gathered in section O 

The vector nonlinear Schrodinger equation (VNLS) on the half line is defined as the 
following initial-boundary value problem for the n-component vector field R{x, t) 



.dR 

dt dx"^ 

R{x,0) = Ro{x) , R{0,t) = go{t) , R,{0,t) = g,{t) , (1.2) 



t— + -2XRR^R = 0, x,tG[0,oo), (1.1) 



where R'^{x, t) is the transpose conjugate of R{x, t) and A is the (real) coupling constant 
which can be normalised to A = ±1. We assume that the functions Rq, Qq and gi live 
in appropriate functional spaces so as to ensure that the calculations in this paper are 
meaningful (for details in the scalar case see for instance the "rigorous considerations" 
section in [18]). In particular, they are decaying functions at infinity. A powerful method 
has been developed by Fokas and various co-workers over the last two decades to tackle 
integrable nonlinear equations on the half-line (see e.g. |12] and references therein). The 
key ingredient is the simultaneous treatment of the two equations appearing in the Lax 
pair formulation of the nonlinear equation. We follow this principle in this paper. 

2 Vector NLS on the line 

In this section, we review the inverse scattering method for VNLS on the line. We follow 
chapter four of the textbook [T9] and collect the main results and notations needed in 
the rest of the paper. 

2.1 Lax pair 

Given R{x,t), define Q{x,t) as the following (n + 1) x (n + 1) matrix-valued field 

(Afitt,^) 

Equation 01. ip is the compatibility condition of the following linear problem for the 
matrix- valued field ^{x,t, kj^ 

$, + ^A;[S3,$] = g<I>, (2.2) 

$t + 22A;2[E3,<I'] = G<I>, (2.3) 



where 



S3 = ( ) and G = 2kQ - iQ^E^ - zQ^E^ , (2.4) 



^ From now on. we drop the x, t (and k) dependence for conciseness unless there is ambiguity. 
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/„ being the n x n identity matrix. Using standard arguments, one can show that 
det ^{x,t,k) = f{k) for all x, t G M, and all values of A; G C for which $ is defined, 
where / can be determined for instance by fixing the asymptotic behaviour of $ as 
X — )■ ±00. Also, if $ and ^ are two solutions of (12.2112. 3D then 



^(x, t, k) = ^(x, t, k) e-^<^(^'*''=)^3 T{k) e*'^^^'*''^)^^ , (2.5) 
where T depends on k only and (j){x, t, k) = kx + 2k'^t. Finally, from 

WQW'^ = -Q^ , (2.6) 

where 

W = (-^/" , (2.7) 

one obtains that W^"^ {x,t, k*)W~^ and $^^(x, t. A;) satisfy the same equations and are 
thus related as in fl2.5p . 

2.2 Direct scattering problem 

We define two fundamental (or Jost) solutions M± of fl2.2|2.3p satisfying 

lim e^'^(^'*'^)^W±(x,t,A;)e-^'^(^'*''=)^3 = 4+1 , keR. (2.8) 

They have the following Volterra integral equation representations 

M±(x, t, k) = 4+1 + r e-*'=(^-^)^«Q(i/, t)M±(|/, t, k)e'^i-y)^^ dy . (2.9) 

J±oo 

Their analyticity and boundedness properties then follow as in the scalar case. Splitting 
the matrices into left and right "column vectors" [§, one can indicate the domain of the 
complex k plane where they are bounded and analytic using the following notations 

M+(x, t, k) = (N{x, t, k), N{x, t, k)) , M_(x, t, k) = (M(x, t, k),M{x, t, k)) , (2.10) 

where A^, M (resp. A^, M) are bounded analytic in the lower (resp. upper) half plane 
as functions of k. From (12. Sp . one finds that 

detM±(x,t,A;) = 1. (2.11) 

Having defined the fundamental solutions, one can proceed with the definition of the 
scattering data. Recalling (12. 5p . we define the matrix S{k) for G M as 

e^'^(^'*''=)^W_(x,t, k)e-'^(^'''''^^^ = e^^(^'*'^)^W+(x,t, k)e-''''^'''''''^^' S{k) . (2.12) 

Several important properties follow from this definition. First, det S{k) = 1. Then, S{k) 
has the following integral representation 

POO 

5(A;) = 4+i+ / e^'y'''-^Qiy,0)M4y,0,k)e-"'y^'dy. (2.13) 



^Here, the left column vector is made of the first n left columns and the right one is made of the 
remaining column. 
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Splitting S{k) into block matrices of natural size^ 

^(^■) - ( ) ■ i^-") 

one gets that a{k) can be analytically continued to Im/c > while a{k) can be ana- 
lytically continued to Imk < 0. It is known that a{k*) = {det a{k))*. In the focusing 
case (A = —1), the usual assumption is that deta(A;) has a finite number of simple 
zeros located in the upper half-plane off the real axis. Denote them kj, Imkj > 0, 
j = 1, . . . , J. Consequently, a{k) has the same number J of simple zeros k* , j = 1, . . . , J 
in the lower half-plane. These zeros play an important role in the inverse scattering 
problem described below. In particular, the residues of b{k)a~^{k) ai k = kj, called 
norming constants, are used to reconstruct the solitonic part of the vector field R. 

2.3 Inverse scattering problem 

Let us rewrite the relations (12.121) as 

M{x,t,k)a-\k) = N{x,t,k)+e^"''^'''''''^N{x,t,k)p{k), (2.15) 
M{x,t,k)d'\k) = e-^'^^'''''''W{x,t,k)p{k) + N{x,t,k), (2.16) 

where p{k) = b{k)a^^{k) and p{k) = h{k)a~^{k). Next, apply the Cauchy projector P~ 
(resp. P+) to the first (resp. second) equality where 

P±(/)(A;) = ^ r ljtl_^du. (2.17) 

Using standard contour integration techniques and the asymptotic behaviour in k that 
can be obtained using integration by parts on the integral representations (12. 9p and 
( I2.13p . one obtains 



o) + ^ iV(x,t,fc,)fi 



^ r e---^---Ar(.,t,.)p(.) 

2tnJ_^ u-ik-tO) ^ ^ 



/Q\ -2ik*x-4ikft_ _ 

N{x,t,k)=(Jj + J2 Nix,t,k*)C, 



I e-^''^^-'^''^^^N{x,t,uj)p{u) 
2m u-{k + zO) 



duj, (2.19) 



where 



Cj = Res{b{k)a-\k), kj) and Cj = Res{b{k)a-\k), k*) . (2.20) 



Note that (12. 6p implies that Cj = ACj. Evaluating the first equation at k = k* and the 
second at A; = kj, one obtains a closed system of linear algebraic integral equations that 



•^For instance, a is n x n and a is a scalar. 
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can be solved in principl^. Below, we will see that this system of equations allows for 
closed-form formulas in the purely solitonic case i.e. when one sets p{k) = = p{k) for 

In the scalar case, the norming constants 7^ are defined as the proportionality con- 
stants between the two column vector forming a Jost solution when the spectral param- 
eter equals to a zero ki 

M{x, t, ke) = e2^<^(^'*''=^)A^(x, t, k^h^ . (2.21) 

When working in the class of exponentially fast decreasing functions, there is a nice 
characterization of 7^ in terms of b{k) which can be analytically continued up to the 
strip {fc G C; < Im A; < K} where K controls the decrease of the functions. With 
K > max{Im kj;j = 1, . . . , J}, we then have 7^ = b{k£). In the vector case, the norming 
constants can be defined as the vectors describing the linear dependence of the column 
of the Jost solution at the zeros kj 

M{x, t, kj)rj = e^"^^'''''''^^N{x, t, kj)Cj , (2.22) 

where Tj is an x matrix and Cj a row vector. Again, in the class of exponentially 
fast decreasing functions one has, without loss of generality, 

r, = A{k,y and Q = -J^6(A:,)A(A:,)* , (2.23) 

where a {k) = '^^'^^^^''^ and A(ky is the transpose of the cofactor matrix of a(k). Al- 
though not the most general, we will work with this convenient characterization to derive 
in particular proposition (13. 2p below. We stress that the results would still be valid in a 
more general case, but this is enough for our purposes of studying soliton solutions. 
Once and A^ are known, one obtains R{x,t) by the formula 

R{x, t) = 2iX J2 e-^^'^^-^^'^f *iV^"^\x, t, k*)C] + - T e-2--4-^*iV^"P\x, t, uj)p{uj) du , 

^ J —00 

(2.24) 



where A^*""^"* is the n x n upper block of A^. This is obtained by comparing the 0(1/ k) 
term in the expansion obtained from (12. 9p by integration by parts and the 0(1/ k) term 
obtained from (I2.19p . 

In the pure soliton case, (I2.18|2.19p yield the following coupled equations 

N{x, t, K) - A i:7V^(^' ^' kl)ClC, . = In , i=l,...,J .(2.25) 

Introducing the matrix p{x,t) whose entries are the following n x n matrices 

Pmi{x, t) = 6meIn-XYl -JT* — h\(h tTT*^-^^' ' (2-26) 



^Tliat this is the case is usuahy estabhshed by using the equivalent formulation in terms of a Riemann- 
Hilbert problem. 
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and the notation N{x,t, kf) = Ni, equations (I2.25P can be compactly rewritten as 

(iVi N2... Nj) t) ={lnln... In) • (2.27) 

One then obtains 

Rix, t) = 2iX {In In... In) ^'\x, t) 



(2.28) 



This is a central formula to understand the behaviour of solitons. In the next section, 
we show how to modify this formula to take into account the presence of the boundary. 
This is then used in Section H] for numerical simulations. 

3 Vector NLS on the half-line 

Although the most general method for tackling initial boundary value problem exactly is 
the one developed by Fokas and co-workers, since we restrict our attention to a family of 
integrable boundary conditions (the vector generalisation of the Robin boundary condi- 
tions), we employ the mirror image technique as used in [17j. These two approaches give 
equivalent results for this class of boundary conditions which would be called lineariz- 
able in the terminology used by Fokas. The advantage of the mirror image technique 
is that it allows one to use the results on the line with minor (but important) modi- 
fications to get explicit solutions. One inconvenience though is that one has to prove 
that the boundary conditions are indeed satisfied by the field. This is because one uses 
the framework presented in the previous section for a modified field P{x,t,k) instead 
of Q{x,t) whose structure is chosen a priori. We note that, to the authors' knowl- 
edge, such a check has only been performed explicitly in [l6j and only in the case of 
Dirichlet boundary conditions. Below, we present a similar proof for the mixed Neu- 
mann/Dirichlet boundary conditions in the vector case. We also show that the class of 
boundary conditions we consider is integrable in the sense that it allows for an infinite 
number of conserved quantities. This is done by exhibiting a generating function for 
such conserved quantities. 



3.1 Deriving integrable boundary conditions 

In order to start with an appropriate form for P{x, t, k), we use the Backlund transforma- 
tion method presented in [13] and further developed in f Til [T5] for the scalar nonlinear 
Schrodinger equation to treat a class of integrable boundary conditions. Contrary to 
the scalar case, in the vector case it is important to study both the x-part and the 
t-part of the equations defining the Backlund transformation to completely characterize 
integrable boundary conditions. This is done below. This is all the more important as 
the t-part of the Backlund method provides the link between this method and Fokas' 
approach in the case of linearizable boundary conditions. 

We introduce a Backlund matrix L{x, t, k) relating the auxiliary problem (12. 2112. 31) for 
$ to the same auxiliary problem for $, with the potential Q replaced by a new potential 
Q, by the equation 

^{x,t,k) = L{x,t,k)^{x,t,k) . (3.1) 
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It is well-known that L, also known as a gauge transformation of the auxiliary problem, 
satisfies the following equations 

L, + tk[^3,L] = QL-LQ, (3.2) 
Lt + 2ik'^[E3,L]=GL-LG, (3.3) 

and preserves the compatibility condition (or zero curvature condition). We look for a 
solution in the form 

L{x,t,k) = kln+i + A{x,t) , (3.4) 
under the symmetry constraint Q{x,t) = Q{—x,t). Writing A in natural block form 



Ai{x,t) A2{x,t) 
A3{x,t) A4{x,t) 



(3.5) 



and inserting in (13. 2p yields the equations 

2^^2(0;, t) = R{-x, t) - R{x, t) , (3.6a) 
-2iA3{x, t) = \ {R\-x, t) - R\x, t)) , (3.6b) 



and 



Ai^.(x, t) = R{-x, t)A3{x, t) - \A2{x, t)R\x, t) , (3.7a) 

A2xix, t) = R{-x, t)A^{x, i) - Ai{x, t)R{x, t) , (3.7b) 

A3^{x,t) = \[R\-x,t)Ai{x,t) - Ai{x,t)R\x,t)\ , (3.7c) 

A4^(x, t) = XR^{-x, t)A2ix, t) - A3(x, t)R{x, t) . (3.7d) 



Note that 



A3{x,t) = \4ix,t) . (3.8) 
Combining (13. 6p and ( ]3.7b|l3.7cp at x = yields the following boundary conditions 

R^{0, t)=- i{A^{0, t)h - Ai(0, t))i?(0, t) , (3.9a) 
i?i(0,t) =ii?t(0,t)(Ai(0,t) -A4(0,t)/„). (3.9b) 

To ensure the compatibility between (I3.9ap and (I3.9bp . we impose 

Ai(0,t) - A4(0,t)/„ = -(Ai(0,t) - A4(0,t)/„)^ , (3.10) 

showing that H = — z(Ai(0, t) — ^4(0, t)/„) is a hermitian matrix. The boundary condi- 
tion now reads 

i?^(0,t) + //i?(0,t) = 0. (3.11) 

We note that at this stage, we have boundary conditions that depend on time a priori. We 
remove this time dependence by requiring that v4i(0, t) and ^4(0, t) be time-independent. 
Then, H is time independent. We see that (13.111) is the vector generalisation of the usual 
Robin boundary condition in the scalar case (r^(0,t) + ar(0,t) = 0, a G M). The fact 
that H is hermitian is the analog of a being real. Let us denote ^4(0) = /3. What we 
have obtained so far reads 

L(0, fc) = A:/„+i + + ° , (3.12) 



with L independent of t at x = 0. To complete the characterization of L(0,/c), insert 
fl312|) in (Q to get 

G(0, k) L(0, k) - L(0, k)G{<d, k) = Q . (3.13) 

Noting that G(0, k) = £3^(0, -A;)S3, this reads 

G(0, -A;) S3L(0, A;) = £3^(0, k)G{0, k) . (3.14) 

It is now convenient to use the hermiticity of H that guarantees that H is diagonahsable 
by a unitary matrix V 

H = VDV\ where D = diag{rfi, . . . , dn}, (3.15) 

and the d/s are all real numbers. Note that the transformation R{x,t) H- V'^RlXjt) = 
R'{x,t) leaves the nonlinear Schrodinger equation invariant and in the new basis the 
boundary condition takes the simple, diagonal form 

R'^{0,t) + DR'{0,t) =0. (3.16) 

This shows that, in the presence of a boundary described by H, the VNLS equation has 
a preferred polarization basis determined by the boundary. In the following, we work in 
this basis and drop the Then, 

L(0, A:) = + ( + ° ) , (3.17) 

and inserting the relation Rx{0,t) + DR{0,t) = in G{0, k), equation (I3.14p yields 

{2i/3In- D)DR{0,t) =0, (3.18a) 
R^{0,t){2i(3In- D)D = 0, (3.18b) 
RR^O, t)D = DRR^{0, t) . (3.18c) 

The compatibility of the first two equations imposes that /3 be purely imaginary: /3 = ia, 
a G M, unless -R(0, t) = (this Dirichlet boundary condition for all the components is 
formally obtained when all the dj are infinite). Then, the first equation shows that either 
dj = or dj = —2a or Rj{0,t) = 0. Finally, the last equation reads 

djRjRl{0,t) = dkRjRliO,t) , j,k = l,...,n. (3.19) 

In general, this means that dj = dk = d i.e. D = din is proportional to the identity 
matrix. The particular case Rj{0, t) = for some j requires some attention. In this case, 
either Rjx{0,t) is also zero and dj = d as before, or in general Rjx{0,t) 7^ 0, meaning 
that dj = 00 is different from the common value d and we must have dj + 2a = 0. So this 
case occurs when formally a = —00. To summarize the results, we have the following 
possible integrable boundary conditions 

Rx{0,t) -2aR{0,t) = , aGM, (3.20) 

or a mixture or Neumann and Dirichlet boundary conditions 

R,{0,t)=0 , jeS , (3.21) 
i?fc,(0,t)=0 , ke{l,...,n}\S, (3.22) 
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where is a subset of {1, . . . , ri}. In the rest of this paper, for convenience we will 
assume that 5* is a proper subset of {1, . . . ,n}. This avoids redundancy for the pure 
Dirichlet or Neumann condition which can be recovered from the Robin condition in the 
limits a — 7- oo or a — 7- respectively. 

In terms of L(0, k) the previous results are more conveniently written by considering 

1 



k + ia 

Note that C is completely equivalent to L since a Backlund matrix is always defined up 
to a function of k but it has the advantage of accommodating the a = — oo case. Then, 
the previous two cases correspond to 

fa. 



C{0,k) 



k—ia 
k+ia 



or £(0,fc) 



(3.23) 



where aj = —1, j E S and aj = 1, j G {1, . . . , n} \ 5. 

We note that in [TU], the previous mixed Neumann/Dirichlet boundary condition 
was not found. In section HJ we will see that this is the most interesting one as far as 
interactions between vector solitons and boundary are concerned. 

It is now worthwhile to comment on relation (I3.14p . This is precisely the relation that 
is imposed in the Fokas' method to obtain the so-called linearizable boundary conditions 
(with the identification J^^L^Ojk) = N{k), see e.g. |12)). Thus, for each solution N{k) 
of equation (I3.14p that can be identified as the value at x = of a Backlund matrix 
L[x,t,k), we see that the so-called linearizable boundary conditions are integrable. In 
fact this is shown in section [3l2] where a generating function for the conserved quantities 
is constructed for the class of boundary conditions described by such an N{k). 



3.2 Integrability 



The fact that our boundary conditions derive from a Backlund matrix ensures the exis- 
tence of an infinite number of conserved quantities. This is what we mean by integrability 
in this context. Adapting a well-known argument for deriving a generating function of 
conserved quantities from the Lax pair to our situation, we show this explicitely below. 
Note that we expect integrability in the Liouville sense to hold as well once one equips 
the present system with a Hamiltonian structure. Indeed, Backlund transformations are 
known to have nice canonical properties (see e.g. [20j). 

Proposition 3.1 A generating function for the conserved quantities of the VNLS with 
integrable boundary conditions is given by 

i{I{k)-l\k*)), (3.24) 



where 



tr 



[R{x, t)(r(x, t, k) — r(x, t, —k))] dx 



and r(x,t, k) satisfies the following Ricatti equation 

= 2ikT + - TRY . 



(3.25) 



(3.26) 
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Proof: The proof follows the one presented in which we adapt to the present vector 
case with a boundary. Recall that we consider two copies of the auxiliary problem 
f l2.2f2.3p related by f l3.ip . All the matrices involved in the argument, like $, G or L are 
split in natural blocks 



Zn Zi2 
Z21 Z22 



Z = ^,G,L. 



(3.27) 



where Zn is an n x n matrix, Z12 is an n-component column vector, etc. Define 
r(x,t,fc) = $2i*n and f{x,t,k) = ^21^11 ix,t, k). Then, ([22]) yields flX^ 

for r and the same equation for F with R replaced by R. Also, we can write 



$ii,$ri' = i?r , <i>m$n = Cii + Gi2r 



-1 



(3.28) 



so using (ln<l>)a;( = (ln$)(^ we get (-Rr)^ = (Gn + Gi2T)x where Gij are the appropriate 
blocks of G. A similar relation holds for F. Finally, from (13. ip we have 

^iit = pn + Li2F)<l>n], 

= [{Ln + LuT)t + {Ln + L,2T){Gn + GuT)]{Ln + LuT)-'^n, (3.29) 

which we compare to $114 = {Gn + Gi2F)<l>ii to get 

tr(Gn + GuT) = tr(Gii + G12F) + tr ln(Ln + Lur)t ■ (3.30) 
Swapping the roles of $ and $ and introducing L = we can also obtain 

tr(Gii + G12F) = tr(Gii + ^12?) + tr ln(Ln + Li2f)t • (3.31) 

This allows us to obtain a more symmetric form of the final result. Putting everything 
together, we get the general result 

dfti 



RT{x,t,k)dx+ / RT{x,t^k)dx 
o Jo 

In (Ln(0, k) + Li2(0, A;)F(0, t, k)) - In ( Ln(0, k) + Li2(0, A;)F(0, t, k) ] .(3.32) 



;atr 



Now, for the problem with boundary, the thing to note is that under the reduction 
R{x,t) = R{—x,t), we have 



r{x,t, k) = —r{—x,t, —k) and L{x,t,k) 



k^ + a^ 



S3L(-x,t,-A;)S3. (3.33) 



Finally, at x = 0, L(0, A;) does not depend on t for the class of boundary conditions we 
have derived and ^12(0, /c) = so the right-hand-side in ( I3.32p vanishes. Therefore, we 
have shown 



dtl{k)=0. 

The special form f l3.24p is used to get real conserved quantities. 



(3.34) 



10 



In practice, the conserved quantities are determined recursively by inserting the 
following expansion 

n=l ^ ' 

into the Ricatti equation to obtain 

n-l 

Ti = -XR^ , r„+i = Tnx + TkRTn-k , n > 1 . (3.36) 

fc=i 

As expected from the presence of the boundary, we see that the conserved quantities 
corresponding to even powers of 1/A; do not appear. In particular, the momentum (order 
l/Zc^) is not conserved. 

The connection between linearizable boundary conditions and Backlund transfor- 
mations was noted in [22|. Here, in view of the previous result, this connection al- 
lows one to identify explicitely the infinite set of conserved quantities once a solution 
£3/7(0, k) = N{k) of the "symmetry" relation fl3.14p is known. Of course, this connection 
goes beyond the specific example of VNLS and holds for any integrable system charac- 
terized by a Lax pair (with possible complications for equations that are not invariant 
under x — )■ —x such as the KdV equation). 

3.3 Mirror image construction 

Following |16[ [T7] . we can now use the result of the previous section to introduce the 
following extended potential 

P{x, t, k) = e{x)Q{x, t) + e{-x)B{k)Q{-x, t)B{-k) , (3.37) 

where 6 is the Heaviside function and 

S(fc)^(^^^) _°i)=S3£(0,fc) 

Note that 

B\k*)B{k) = In = B{-k)B{k) , keC 

We have the property 

P{-x, t, -k) = -B{-k)P{x, t, k)B{k) . (3.40) 

The idea is to apply the inverse scattering method as explained in section [2]to P instead 
of Q. By construction, for x > 0, any solution found for P gives a solution Q of VNLS 
on the half-line. The delicate point is to check whether this solution does indeed satisfy 
the boundary conditions that we derived in the previous section and that we encoded in 
Fourier space in the matrix B{k) obtained from £(0, k). 

For conciseness, we keep the same notations as in the previous section but with Q 
replaced by P. Using iK^ . we obtain that M+{x,t,k) and B{k)M_{-x,t, -k)B{-k) 
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(3.35) 



(3.38) 



(3.39) 



satisfy the same Volterra integral equation. Noting that they have the same asymptotics 
as a; — 7- oo we deduce 

M+{x,t,k) = B{k)M4-x,t,-k)B{-k). (3.41) 

From this and fl2.12p . we obtain 

S-\k) = B{k)S{-k)B{-k) . (3.42) 

Next, from ([22]), we have that WMi{x,t, k*)]^^ = M^\x,t, k) which in turn yieldfl 

S'\k) = WS\k*)W-\ (3.43) 

Combining the two previous equations yields the important result 

WS\k*)W'^ = B{k)S{-k)B{-k) . (3.44) 

In components, this reads 

a^{k*) = B{k)a{-k)B{-k) , d*{k*) = a{-k) , (3.45) 
b\k*) = \h{-k)B{-k) , \h\k*) = B{k)b{-k) . (3.46) 

Note that due to the property B\k*) = B{—k), the last two relations are consistent. 
The first relation implies that if kj is a zero of deta(A;) then —k* is also a zero. The 
same holds true for a{k). Therefore, we find that the main observation of [17\ that the 
relevant zeros to formulate the inverse problem come in quartets {^kj, ^k*} is also valid 
in the vector case. The total number of zeros in each half-plane is even: J = 2N and 
there are zeros in each quadrant of the complex plane. For each j = 1, . . . ,N, it is 
enough to assume that kj is in the first quadrant i.e. Re kj > and Im kj > 0. The other 
three zeros are obtained by successive reflections with respect to the imaginary and real 
axes. 

Associated to these zeros are norming constants which we denote Cj (for kj) and C/ 
(for —kj). From the symmetry (12.61) . the norming constants associated to k* (resp. —kj) 

are then ACj (resp. AC^-^). 

The subtle point of the method now is to reconstruct R{x, t) from the previous 
scattering data corresponding to P{x,t, k). This is achieved by establishing the analogs 
of (I2:T8|) . (I2:T9|) and flXM]) and restricting to x > where P{x,t,k) = Q{x,t). To 
do so we must investigate the analytic properties of M±{x,t,k) as functions of k for 
X > 0. The argument follows that given in [16]. For x > 0, P{x,t,k) = Q{x,t) 
so the analytic properties of M_|_ are the same as in the infinite line case. Then, the 
analytic properties of M_ are obtained using the symmetry relation (13.411) . In the 
mixed Neumann/Dirichlet case, no extra pole arises from the presence of B{k) (which 
is in fact independent of k). In the Robin case, M gets an extra pole at k = ia if 
a > and M gets an extra pole at k = —ia if a > but are otherwise analytic in 
their respective half-plane. Then, a{k) and d{k) share the same analytic properties as 

^We note that strictly speaking, S{k) is defined first for fc G K. But as we have seen, its entries can 
be analytically continued to various parts of the complex plane. Eq p.43p is then to be understood 
as relations for the entries of 5* valid wherever they make sense. Consequently, the same is true of eq 
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M and M. As a consequence, the fundamental quantities Ma~^ and Ma~^ have at 
worst a removable singularity. Therefore, the analysis of section 12.31 goes through as 
before, taking into account the additional symmetry properties of the scattering data 
discussed above. This is summarized in the following proposition which constitutes the 
main result of this work. 

Proposition 3.2 The reconstruction formulae for Q{x,t), a; > are given by 112. 18\} . 
112. 19\} and 1^2. 24\ ) with J = 2N and, without loss of generality, the substitution ^2^-1 
kj, j = 1, . . . , N with corresponding norming constants Cj and — )■ —k*, j = 1, . . . ,N 
with corresponding norming constants C'j. The scattering data involved in the formula 
satisfy the symmetry relations (3^J^ \3.46^ and 



XC';C, = B{-k,)^^ , j = l,...,N. (3.47) 



Proof: The only thing left to prove are the relations fl3.47p between the norming constants 
Cj of the real solitons (those living on x > 0), and the norming constants Cj of the mirror 
solitons (those living on x < 0). We use the following chain of relations obtained from 
fl2.12p and (13.41 p evaluated at A; = kj and k = —kj 



M{x,t,kj)A\kj) = e^''t'^'''''''^^N{x,t,kj)b{kj)A\kj) , (3.48) 

N{x,t,kj) = 1^ ^^M^-x,t,-kj), (3.49) 

M{-x,t,-kj) = e'^'^'-''''''~^^^N{-x,t,-kj)b{-kj), (3.50) 

N{-x,t,-kj) = 1^ M{x,t,kj)B{kj) . (3.51) 



The compatibility of these relations yields A^{kj) = —B{kj)b{—kj)b{kj)A^{kj) which 
is just fl3.47p upon using Cj = ^,^^^ b{kj)A{kjy, Cj = ^^^^^-^ = ^Cj and a {—kj) = 



-a'{kj). 

Note that following the same argument but for k = —k* and k = k*, we arrive at the 
following similar relations 

ACjc; = 5(/c;)^tS)! , j = l,...,N, (3.52) 

which are compatible with fl3.47p . 

The symmetries of the scattering data ensure that Q satisfies the desired boundary 
conditions. We show this in general in the mixed Neumann/Dirichlet case. 

Proposition 3.3 Without loss of generality take B{k) = diag{l, ...,1,— 1,...,— 1) = 

B then the reconstruction formulae yield a potential R which satisfies 

R{-x,t) = BR{x,t) , (3.53) 
i.e. R satisfies mixed Neumann/Dirichlet boundary conditions. 
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Proof: In the case of fast decreasing solutions, the reconstruction formula fl2.24p can be 
compactly rewritten as 

R{x, t) = - [ e-2^""-^'"'*iV^"''V, t, uj)p{uj) du , (3.54) 
Jc 

where C is a contour in the lower half-plane from — oo to +00 that passes below all the 
zeros of a{k) and that is symmetric with respect to the imaginary axis. Repeating the 
arguments of section 12.31 from the relation 

e^<l>i^,t)^3M_{x, t, fc)e-^*(^'*)^^ S-\k) = e^'^(^'*)^W+(x, t, A;)e-*'^(^'*)^3 , (3.55) 

instead of ( I2.12p . one can derive an equivalent reconstruction formula 

R{x,t) = -- [ e-2'""-^'"'*M("P)(x,t,a;)r(u;)du;, (3.56) 
Jc 

where M (and M) satisfy coupled integral equations similar to fl2.18|2.19l) and T{k) = 
d{k)c~^{k) is obtained from the entries of S^^{k) 

. _ f c{k) d{k) 
^ ^"f-y d{k) c{k) 

Now combining f l3.43p and (I3.44p . we get r(A;) = —Bp{—k). Also, from (I3.4ip . we get 

{—x^ t, —k) = BM^'^P\x, t, k)B. Finally, performing the change of variable u — )■ —u 
in (I3.54P and using the previous relations yields the result. ■ 

To the authors' knowledge, a general proof that the field obtained with the mirror 
image technique obeys the desired boundary conditions is not known within this for- 
malism. Even in the scalar case, an explicit argument has been given in only in 
the case of Dirichlet boundary conditions. Arguments based on the linear limit of the 
nonlinear Schrodinger equation in the case of general boundary conditions have been 
given in [T7]. It seems that the proper way to attack the problem is to use the general 
method developed by Fokas and co-workers in the case of linearizable boundary condi- 
tions. In view of the connection that we have made above between this method and the 
Backlund transformation approach that we have used here in the vector case, we see 
that a full treatment of the problem requires to formulate the vector version of Fokas 
method. Hopefully, this problem is still amenable for the class of linearizable boundary 
conditions that we consider here. This is beyond the scope of the present paper and is 
left for future investigation. 

It is worth noting the appearance of the matrix A^[kj) in the relations between the 
norming constants of the "real" solitons and those of the mirror image solitons. This is a 
new feature compared to the scalar case. Using the dressing method [23], it can be shown 
that these matrices can be expressed in terms of all the norming constants. Hence, the 
system (I3.47P is a system of coupled matrix equations that one has to solve in order the 
completely characterize the soliton solutions with boundary. This is technically more 
challenging than in the scalar case where one could use directly the trace formula for 
a{k) to solve the problem [T7]. This is reminiscent of the difference between the scalar 
NLS and the Manakov model on the line where the presence of polarization vectors for 



(3.57) 
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the solitons greatly complicates the analysis of the soliton interactions [9]. Here, the 
presence of the cofactor matrix captures the influence of the polarization vectors of the 
solitons on their interaction with the boundary. 

A general efficient algorithm has not yet been found to solve these coupled equations 
and completely characterize the A^-soliton solutions with a boundary. We hope to return 
to this question in a future work. In the next section, we show how this can be done for 
one soliton and use the result to investigate its interaction with the boundary. 



4 One soliton reflection: transmission between modes 

We consider the pure soliton case with = 1 z. e. we study the behaviour of one soliton 
bouncing off the boundary. The scattering data is ki and Ci. From this, we deduce the 
scattering data corresponding to the mirror soliton: —k^ and determined by (13.521) 
(with A = -1) 

(4,1) 

We use the dressing method to compute A{—kiy. Let us first recall this method to 
generate two arbitrary solitons from the vacuum solution and then we will apply it in 
the present special case. Let Hi, k,2, Ci and C2 be the poles and norming constants for 
the two solitons. The scattering coefficient a{k) is then given by 

a{k) = D,{k)D2{k) , (4.2) 

where 

D,{k) = I^+(^-l]u,, (4.3) 



and Ilj is a (one-dimensional) orthogonal projector obtained recursively by requiring 
Cj a{Kj) = 0, j = 1, 2. From this, one easily deduces 

= - ^;)U2D,{K2r\ (4.4) 

with CiIIi = Ci and p2^2 = P2 where p2 = C2-Di(/t2)- 

Now, we insert this in (14.11) with the identification Ki = ki, k,2 = —kl, Ci = Ci and 
C2 = C'l- Multiply on the right by Di{k2) and project on p\ to obtain 



lB{-kl)Cl^Vl (4.5) 



p\ ^ a'{-kl 
P2P2 ^1 ~ ^1 
Hence, we deduce the solution 

For practical calculations, we note finally that a'{—kl) can be obtained from the well- 
known trace formula for det a{k) (see e.g. [19j) which gives, in our case. 



C[ = ^D,{-kl)-K (4.6) 



a\-kl) = -h±tL^ _ (4.7) 

^ 2kl{ki-kl) ^ ' 
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In the following, we present numerical results for the case of two components {n = 2) 
and for the mixed Neumann/Dirichlet case. The Robin boundary condition does not 
bring anything new in the vector compared to the scalar case in the sense that 

in any polarization basis, each component of R always satisfies a scalar Robin boundary 
condition. However, in the mixed Neumann/Dirichlet case, if the polarization basis is 
different from the preferred boundary basis i.e. if we restore the unitary V matrix in 
fl3.15p . then an interesting phenomenon of reflection-transmission between the modes 
appears: the amplitude of the soliton envelope \Rj{x,t)\ is different before and after its 
interaction with the boundary! 

As an example, we take 

Y _ f cos 9e^^ sin Oe^^ 
y — sin 6'e~*^ cos 0e~*^ 

with 6', ^, G M and consider the boundary matrix Bg given by 

Be = v(^\ \ ^V~\ (4.9) 

The parameters 6,C,C, measure the "deviation" from the natural boundary basis corre- 
sponding to 6* = C = ^ = 0. Below are plots of \Ri\ and |i?2| in x, t space for /ci = 1 + | 
and Ci = (2 1)*. 

First, for 6 = ( = ^ = (Fig. [1] and [2]) we see that, as expected, the first mode 
satisfies a Neumann boundary condition while the second mode satisfies a Dirichlet 
boundary condition. These plots are very similar to those of [T7] in each case. Each 
mode behaves like a scalar solution seeing its own boundary condition. 





Fig. 1: Reflection of the first component of the 
soHton on the boundary with Neumann bound- 
ary condition. The ampHtude before and after 
the reflection is the same. 



Fig. 2: Reflection of the second component 
of the soUton on the boundary with Dirichlet 
boundary condition. The amplitude before and 
after the reflection is the same. 



Now, for ^ = I and ( = ^ = (Fig. [3]and|4j), we clearly see that the amplitude of 
each mode before and after the interaction with the boundary is different. Both modes 
are reflected but part of mode 1 is transmitted to mode 2. 

This shows that the boundary acts as some sort of polarization filter. We emphasize 
though that there is no loss in the transmission process in the sense that the quantity 
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Fig. 3: Reflection of the first component of the 
sohton on the boundary. The outgoing mode has 
a smaUer amphtude than the incoming one. 



Fig. 4: Reflection of the second component 
of the sohton on the boundary. Part of the 
first mode has been transmitted to the outgo- 
ing mode here. 



+ I-R2P is indeed a conserved density. Another way to see this is to investigate 
the same quantity in the reference frame of the incoming or outgoing sohton. Here the 
velocity of the incoming (outgoing) soliton is —v {v) with f = 4Re/ci = 2. We find 
indeed that 



lim + \R2\^){-vt + Xr 



t— >— 00 



lim(|i?i| 



+ \R2\^){vt + Xout.t) 



(4.10) 



where Xin {xout) is the incoming (outgoing) soliton centre. 

One can then wonder if one can use the boundary to control the amplitudes of the 
modes by changing the values of the boundary parameters. The following plot shows 
that this can be done. We keep fci = 1 + | and Ci = (2 1)* and we vary 6 from to | 
while fixing ( = 1.11 and ^ = 0. The incoming amplitudes are of course constant (lines) 
while the outgoing amplitudes vary (dashed curves). In our example, we see that the 
second mode (in blue) can be made vanishingly small (for 6 ^ 1.15). The total outgoing 
amplitude is constant though (black line) and equals the total incoming amplitude. 
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Fig. 5: Amphtudes in the sohton frame as functions of 0. The red hne is the incoming 
amplitude. The blue line is incoming |-R2|- Dashed curves are the corresponding outgoing 
amplitudes. The black line is the total amplitude Y^|i?ip + |-R2p- 
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5 Conclusions 



The main goal of this paper was to estabhsh the inverse scattering formulae for the 
VNLS model on the half-line with integrable boundary conditions. We had first to 
identify such integrable boundary conditions and this was done using the Backlund 
matrix approach. An important difference with the scalar case is the need to consider 
the time part of the Backlund transformation. Indeed, for the scalar case, this does 
not lead to new information and, in fact, it was not considered this way in the original 
papers [T^ [T5] . Only the space part was considered and then it was checked that the 
solution was compatible with time evolution. Here, we considered the time part directly 
and found that the solution to the space part has to be further reduced to two classes of 
possible integrable boundary conditions (which are then automatically compatible with 
the time evolution). On the technical side, this is chiefly due to the fact that in the 
scalar case, there is no difference between RR^ and R^R while in the vector case, the 
former is a square matrix of size n while the latter is a scalar. 

There seems to be a sharp contrast between vector and scalar cases in terms of which 
boundary conditions holding in the linear case go through to the interacting case. In- 
deed, in the scalar case, the Robin boundary condition qx{0,t) + dq{0,t) = 0, easily 
established for the free Schrodinger equation survives the switching on of the interaction 
in the sense that it is compatible with the integrability of the nonlinear Schrodinger 
equation. But here in the vector case, only subfamilies of the general Robin condition 
Rjx{0, t) + djRj{0, t) = 0, j = 1, . . . ,n survive the requirement of integrability of the non- 
linear equation when implemented in the form of the Backlund transformation method. 
Interestingly, this is not seen if one only looks at the x-part of the Backlund transfor- 
mation method and the analysis of t-part also yields an important connection with the 
notion of linearizable boundary conditions in the Fokas method. In the Appendix, we 
examine in more details the linear and nonlinear systems. 

Using the boundary conditions thus found, we implemented the vector version of the 
mirror image technique [T7] to obtain the reconstruction formulae for the fields. Here 
again, a major difference with the scalar case is the appearance of coupled equations to 
determine the norming constants of the "virtual" solitons from the given scattering data 
characterizing the "real" solitons. We solved this problem explicitely in the case of one 
soliton. Using these results, we discovered an interesting phenomenon: the boundary can 
be tuned to alter the polarization of a soliton. This was demonstrated with numerical 
simulations. 

The question of an efficient treatment of the general N-soliton problem in the presence 
of a boundary is an exciting open issue. 
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Appendix: Similarities and differences between linear 
and nonlinear Schrodinger equations with a boundary 

We take some time to discuss the classes of integrable boundary conditions that we 
have found using the Backlund transformation method. It is well known that for the 
free vector Schrodinger equation on the half-line (corresponding to switching off the 
nonlinearity by setting A = 0), the following general boundary conditions hold 

{U - In)R{0, t) + l{U + In)RM ^) = , (5.1) 

where U is some n x n unitary matrix. This is the same as (13.111) (and hence (I3.16P 
in the boundary basis) where we allow formally some of the eigenvalues of H to be 
infinite if the corresponding eigenvalues of f/ + /„ are zero. But with the nonlinear 
term, we found that only subclasses of these general boundary conditions survive if 
we want to maintain integrability. Therefore, it appears that the nonlinearity imposes 
severe restrictions. One can then wonder if such restrictions are due to the use of the 
Backlund transformation method and if they could be circumvented by a generalization 
of it. There is an indication that the restrictions are genuine or at least consistent 
with existing methods. Indeed, in the quantum version of the system considered here, 
the integrable boundary conditions are controlled by the reflection equation [4j. The 
solutions for the corresponding reflection matrix have been classified in [7j and the same 
reduction of possible boundary conditions appears due to the nonlinearity. Algebraically, 
this means that the R matrix of the Yangian of gl{n) appears in the reflection equation 
instead of just the identity matrix in the free case. 

To conclude this discussion, it is instructive to perform an analysis of the free case 
along the line of the Lax pair approach and compare with the nonlinear case. This sheds 
some more light on the differences between the two cases and on why integrability of the 
nonlinear case is likely to impose severe restrictions on the allowed boundary conditions. 
In doing so, we provide an alternative and more direct proof of proposition 13.11 in the 
mixed Neumann/ Dirichlet case. In particular, this proof does not involve the Backlund 
matrix L. 

Keeping the notations of section 12. H a Lax pair formulation for the free vector 
Schrodinger equation 

iQt - g..S3 = , (5.2) 

reads 

+ ik[E3, fx] = Q , (5.3) 
fit + 22A;2[S3, fi] = 2kQ - tQ^J:^ . (5.4) 

Splitting as usual in natural blocks, the linear analog of F in proposition 13. II is /i2i = r° 
which satisfles the following linear version of the Ricatti equation (I3.26P 

r° = 2tkT^ + R^ , (5.5) 



which gives a series expansion r° = . " with r° = —d!^~^R^. We also have 

n=l ^ ' 

(r°i?)t = i(r^Rx — T^R)x and hence, on the half-line, a generating function for the 
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conserved charges is 

POO 

I\k)= {T\x,t,k) -T^{x,t,-k)) R{x,t)dx, (5.6) 
Jo 

with 

dti\k) = i (r°(o, t, k) - r°(o, t, -k)) R{o, t) - i (r°(o, t, k) - r°(o, t, -k)) r,{o, t) . 

(5.7) 

Collecting everything, the charges are 

/°= / dlPR^Rdx , p>0, (5.8) 

and they are conserved iff 

{dlP+^R^R-dlm^d^R)^^^ = Q , p>0. (5.9) 

Looking at p = 0, the usual charge /q is conserved when the current is zero at a; = 0: 
[R^R — R^ Rx) = 0. This is equivalent to (15. ip and this is how these boundary 
conditions arise in the first place (note that they also ensure that the energy is real). 
It is rather easy to adapt the Fourier transform method to incorporate these bound- 
ary conditions. The beauty of the linear case is that then, all the higher charges are 
automatically conserved (this is precisely the main difference with the nonlinear case). 
Indeed, the general solution to the free Schodinger equation with boundary conditions 
(15. ip can be written as 

R^x,t)= —e''''-"'''A{k), (5.10) 

where A is a vector-valued function with the property A{k) = B{—k)A{—k) where 

B{k) = -[{U- In) + k{U + InT' [{U - In) - k{U + /„)] . (5.11) 

But then, we get more than just (15. ip . In fact, the reconstructed field (I5.10p satisfies 
the following mirror image symmetry 

[{U - In)R + l{U + In)R.] (X, t) = - [{U - In)R + l{U + /„,)i?x] ("X, t) , (5.12) 

of which ( 15. ip is a consequence. Hence, 

[{U - Qdl^R + i{U + In)d^J'^^R] (x, t) = - [{U - Qd^^R + i{U + Qd^'^^R] (-x, t) , 

(5.13) 

which immediately implies (15.91) as required. 

Now for the nonlinear case, using the ingredients in the proof of proposition 13. H we 
have 

dt {T{k)-V{-k))Rdx = [i{Tx{k)-V,{-k))R-i{V{k)-V{-k))Rx 
Jo 

+ i(T{k)RT{k)R - T{-k)RT{-k)R)]^^Q , 

(5.14) 
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with 



= 2ikT + - TRY . 



(5.15) 



Comparing with (15. 7p and ( 15. 5p . the nonhnearity shows up at two levels: first there is an 
extra term in the right hand side of (15.141) and second, F satisfies a nonlinear equation 
with the extra TRT term compared to r°. Despite this fact, an explicit calculation shows 
that for the first charge, everything is the same as in the linear case. So one is led to 
impose (15. ip . A crucial difference now is that one cannot use the Fourier transform but 
has to use the inverse scattering method instead. It is quite plain that this is likely to 
impose further restrictions on the possible integrable boundary conditions. The proof of 
proposition 13.11 shows that the Backlund transformation method selects such integrable 
boundary conditions. Although it is not a general proof, a careful inspection of the first 
few higher conserved charges (already for n = 2) shows that there is little hope that 
more general boundary conditions than the ones we found ensure the conservation of 
all higher charges. In contrast, we now show that for the mixed Neumann/Dirichlet 
conditions, a nice argument similar to that of the linear case can be built to conclude 
that the left-hand side of (I5.14p vanishes. This provides a direct check of proposition 
13.11 and is based on the result of proposition 13.31 i.e. R{—x) = BR{x). The latter fact 
implies the important relation mirror image relation 



which is in fact the nonlinear analog of i?^(— x) = R\x)B together will all the higher 
derivatives of this relation. As a consequence, r(0, t, A;)i?(0, t) = — r(0, t, — A;)-R(0, t), 
r(0, t, A;)i?^(0, t) = r(0, t, -k)R^{0, t) and r^(0, t, k)R{0, t) = 1^(0, t, -k)R{0, t). Insert- 
ing in the right-hand side of (I5.14p yields zero as required. 

References 

[1] S. Manakov, On the theory of two-dimensional stationary self-focusing of electro- 
magnetic waves, Sov. Phys. JETP 38, (1974), 248. 

[2] E. H. Lieb, W. Liniger, Exact Analysis of an Interacting Bose Gas. I. The General 
Solution and the ground state, Phys. Rev. 130 No. 4 (1963) 1605. 

[3] C. N. Yang, Some exact results for the many-body problem in one dimension with 
repulsive delta-function interaction, Phys. Rev. Lett. 19 (1967), 1312. 

[4] E. K. Sklyanin, Boundary conditions for integrable quantum systems, J. Phys. A21 
(1988), 2375. 

[5] S. Murakami, M. Wadati, Gonnection between Yangian symmetry and the quantum 
inverse scattering method, J. Phys. A29 (1996), 7903. 

[6] M. Mintchev, E. Ragoucy, P. Sorba, Ph. Zaugg, Yangian symmetry in the Non 
Linear Schrddinger hierarchy, J. Phys. A32 (1999), 5885. 

[7] M. Mintchev, E. Ragoucy, P. Sorba, Spontaneous symmetry breaking in the gl(N)- 
NLS hierarchy on the half line, J. Phys. A 34 (2001), 8345. 



T{x, t, k) 



V{-x,t,-k)B 



(5.16) 



21 



[8] A. P. Fordy, P. P. Kulish, Nonlinear Schrddinger equations and simple Lie algebras, 
Comm. Math. Phys. 89, (1983), 427. 

[9] T. Tsuchida, N-soliton collision in the Manakov model, Prog. Theor. Phys. Ill, 
(2004), 151. 

[10] I. T. Habibulhn, S. I. Svinolupov, Integrable boundary value problems for the mul- 
ticomponent Schrddinger equations, Physica D 87, (1995), 134. 

[11] M. Ablowitz, H. Segur, The inverse scattering transform: Semi-infinite interval , J. 
Math. Phys. 16, (1975), 1054. 

[12] A. S. Fokas, A Unified Approach to Boundary Value Problems, CBMS-SIAM (2008). 

[13] I. T. Habibullin, Bdcklund transformation and integrable boundary-initial value 
problems, in Nonhnear world", vol. 1 (Kiev, 1989), 130-138, World Sci. Publ., River 
Edge, NJ. 1990; Integrable initial-boundary value problems, Theor. and Math. Phys. 
86 (1991),28. 

[14] P. N. Bibikov, V. O. Tarasov, A boundary-value problem for the nonlinear 
Schrddinger equation, Theor. Math. Phys. 79 (1989), 334. 

[15] R. F. Bikbaev, V. O. Tarasov, Initial-boundary value problem for the nonlinear 
Schrddinger equation, J. Phys. A24 (1991) 2507. 

[16] A. S. Fokas, An initial-boundary value problem for the nonlinear Schodinger equa- 
tion, Physica D 35, (1989), 167. 

[17] G. Biondini, G. Hwang, Solitons, boundary value problems and a nonlinear method 
of images, J. Phys. A42, (2009), 205207. 

[18] A. S. Fokas, A. Its, The linearization of the initial-boundary value problem of the 
nonlinear Schrddinger equation, SIAM J. Math. Anal. 27, (1996), 738. 

[19] M. J. Ablowitz, B. Prinari, A. D. Trubatch, Discrete and continuous nonlinear 
Schrddinger systems, London Mathematical Society Lecture Notes Series 302, Cam- 
bridge University Press, 2003. 

[20] E. K. Sklyanin, Canonicity of Bdcklund transformation: r-matrix approach. II, 
solv-int/9903017. 

[21] V. Caudrelier, On a systematic approach to defects in classical integrable field the- 
ories, IJGMMP 5, (2008), 1085. 

[22] A. S. Fokas, Linearizable initial boundary value problems for the sine-Gordon equa- 
tion on the half-line, Nonlinearity 17, (2004), 1521. 

[23] V. E. Zakharov, A. B. Shabat, Integration of nonlinear equations of mathematical 
physics by the method of inverse scattering. II, Funct. Anal. Appl. 13, (1979), 166. 



22 



